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Abstract 

Using the framework of metriplectic systems on R 71 we will describe a constructive geometric 
method to add a dissipation term to a Hamilton-Poisson system such that any solution starting in 
a neighborhood of a nonlinear stable equilibrium converges towards a certain invariant set. The 
dissipation term depends only on the Hamiltonian function and the Casimir functions. 
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1 Introduction 

In an attempt for an unification of the conservative and nonconservative (or dissipative) dynamics, 
A.N. Kaufman 4J has introduced the notion of metriplectic system. More exactly, let {., .} be a Poisson 
structure on R™ with II the associated matrix i.e., II = [{xi, Xj}] and C\, Ck € C°°(R™, R) a complete 
set of functionally independent Casimir functions. Let G be a smooth function from R n to the vector 
space of symmetric matrices of type n x n. 

Definition 1.1. (J^j) A metriplectic system on R" is a system of differential equations of the following 
type: 

x = U(x) ■ VH(x) + G{x) ■ Vip(C u C k ){x) (1.1) 
where H e C°° (R™ , R) and <p € C°° (R fc , R) such that the following conditions hold: 
(Ml) IT • VCi = 0, i = l,k, i.e. Ci is a Casimir of our Poisson configuration (R", {•, •}). 
(M2) G VH = 0. 

(MS) (Vp(Ci, . . . , C k ))* ■ G ■ Vp(Ci, . . • , C k ) < 0. 
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Remark 1.1. It is not hard to see that: 

dH 

(i) The derivative of H along the solutions of (jl.ip is —j— ~ 0, i.e. H is a conserved quantity of the 
dynamics (jl.ip . 

(ii) The derivative of <f(Ci, . . . , Cfc) along the solutions of (jl.ip is — - < 0, i.e. ip(C±, . . . , Cfc) 

plays the role of the "entropy function" for the dynamics (jl.ip . 

□ 

Remark 1.2. The dynamics (jl.ip can be viewed as a perturbation of the Hamilton-Poisson system: 

x = {x, H}, 

with the dissipative term G ■ V</?(Ci, . . . , Cfc). □ 

Metriplectic systems have been extensively studied in connection with mathematical physics prob- 
lems, see for instance [2], [4], [5], [6] and [8]. In [I] was proven that a certain type of dissipation induces 
instability. The goal of our paper is to introduces a new type of dissipation, using the formalism of 
metriplectic systems, in such a way that any solution of (jl.ip starting in a small neighborhood of a 
nonlinear stable equilibrium converges towards a certain invariant set containing the equilibrium. The 
dissipative part can be interpreted as a set of controls added to the conservative part. 



2 A class of metriplectic systems on MJ 1 

Let (R", {■,■},#) be an Hamilton-Poisson system on 1" and Ci,...,C k £ C 00 ^",]^) a complete set 
of functionally independent Casimir functions of the Poisson vector space (M n , {•,•})• The dynamic is 
described by the following set of differential equations: 

x = n-VH. (2.1) 

Our goal is to construct explicitly a dissipative perturbation, i.e., to determine effectively a matrix 
g = [g lJ ] such the perturbed system: 

± = n- V-ff + G- V¥?(Ci,...,C fe ) 

to be a metriplectic system. 
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def 8H def dC '■ 

Denote with diH = — — and diCj — The matrix G = [g y l given below satisfies all the 

oxi axi 



conditions from the definition of a metriplectic system 

n 

- E (cW) 2 diifd 2 if 



G 



» ^ i 

d 1 Hd 2 H 



E (^ff) 2 



d!Hd n H 



d 2 Hd n H 



i = 1 
» 7 s 2 



E (5^) 2 



(2.2) 



Indeed, we have that the j-component of the vector field G • VH is given by 

n n 

(G-VH), = ]>>^) 2 ^ff + (-]>>^mff 
= 0. 

Consequently, condition (M2) is satisfied. 

In the case when in the dissipation term we take only one Casimir function we have the following 
computation, 



n n 



(VC)< • G ■ VC = J2 d i C ( G ■ VC h = J2 djiJ^idiHdjHdtC - {d^fdjC)) 

= d 1 Cd 2 H(d 2 Cd 1 H - + • • • + <9iC<9„iJ(<9„C<9iff - <9iC<9„ff) 

+ d-iCdxHidxCdrf - d 2 CdxH) + ■■■ + d 2 cd n H(d n cd 2 H - d 2 Cd n H) 

+ : 

+ dnCfrHfrCdnH - OnCdiH) + ■■■ + d n Cd n -lH \d n -lCd n H - OnCdn-iH). 

Regrouping the terms we obtain the desired inequality which proves that condition (M3) is also satisfied, 
(VC)*-G-VC - -{dxCd-iH -8-iCdxH) 2 -{dxCd^H - d^CdxH) 2 -{d x Cd n H - dnCdxH) 2 

-(d 2 cd 3 H -d 3 cd 2 Hf -{d 2 cd n H -d n cd 2 H) 2 



< 0. 



(2.3) 
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Remark 2.1. The above inequality is an equality iff dH AdG — 0. Consequently, we obtain equality in 
(12. 3[) iff ViJ and VG are linearly dependent. If, instead of one Casimir function, we take a combination 
G := if (Ci, . . . , Cfc) of the complete set of functionally independent Casimir functions Ci, . . . , C k we 
obtain a new Casimir function G and consequently, 

(V^Ci, . . . , C fe ))* • G • Vip{Cx, ...,C k ) = (VC) 1 ■ G ■ VG < (2.4) 

with equality iff V-ff and V<p(Gi, . . . , G/c) are linearly dependent. □ 

Let us consider now the metriplectic system (jl.ip where G is given by the relation (|2.2[> . Then we 
can define in a canonical way two vector fields on R™, namely: 

= n • VH 

and 

e = n- vff + G- v^(Ci,...,c fc ). 

Proposition 2.1. Let (K n ,n, i/) fee an Hamilton-Poisson system. If xo G R" is an equilibrium state 
of the vector field £, i.e., £(a;o) = i/ien is an equilibrium state of the vector field £n- 

Proof. Indeed, £,(xo) = implies that 

(Vp(C 1 ,...,C k )(x o )) t t(x o ) = 0, 

and so 

(V^Gr, . . . , C k ){x )Yll(xo)VH{xo) + (V^(C 1; . . . , G fc )(x )) t G(a ; o)V^(G 1 , . . . , G fe )(x ) = 0. 

This is equivalent (since V</?(Gi, . . . , C k ) is a Casimir we have (V<p(Gi, . . . , Gfe)(a;o)) t II(xo)Vif(a;o) = 0) 
with 

(V^Ci, • ■ ■ , G fe )(ar )) t G(x )V^(G ll . . . , C fc )(a; ) = 0. 
This leads us immediately via Remark 12.11 to 

VH(x ) = \V<p(C 1 ,...,C k )(x ) 

for some A el". Therefore 

£n(zo) = U(x )VH(x ) 

= \Tl(x )V<p(C 1 ,...,C k ){xo) 
= 

as required. □ 



Proposition 2.2. Let (K n ,II, H) be an Hamilton-Poisson system. Let xq 6 R™ be an equilibrium point 
of the vector field £n- If there exists a function ip G C°°(R fc ,R) such that V<p(Ci, . . . , Ck){xo) and 
VH(xq) are linear dependent then xq is an equilibrium point of the vector field 



Proof. Indeed, if 



e = n-vir + GV¥»(Ci,...,c fc ). 



then we have also that 

n(» ) • Vfr(ar ) = 
and consequently we have two possibilities: 

(i) WH(xo) = 0. This implies (see the construction of G) the equality 

G{x ) = 

and then 
Therefore 

£(zo) = 0. 

(ii) H(xq) ■ VH(xq) — and VH(xq) / 0. By hypothesis we obtain 

G(x )-\7<p(C 1 ,...,C k )(x ) = XG{x )VH{x ) 

= 

and we can conclude that 

f(zo) = 

as required. □ 

Corollary 2.3. The set E := {x £ R"|V_ff(x) and \7ip(Ci, . . . , Ck)(x) are linearly dependent} is a set 
of equilibrium points for both vector fields £n and £. 

Proof. For an arbitrary point y £ E we have 

n(y) • VH(y) = AIT(y) • V<p(C u C k ){y) = 0, 

which shows that y is a equilibrium point for the vector field £n- Proposition 12.21 implies that y is also 
an equilibrium point for the vector field £. □ 
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3 Metriplectic systems and asymptotic stability 

For the beginning let us briefly recall some definitions of stability for a dynamical system on K™ 

± = /(*), (3.1) 

where / E C°°(R™,R n ). 

Definition 3.1. An equilibrium point x e is stable if for any small neighborhood U of x e there is a 
neighborhood V of x e , V C U such that if initially x(0) is in V , then x(t) E U for all t > 0. If in 
addition we have: 

lim x(t) = x e 

t — >oo 

then x e is called asymptotically stable. 

For studying more complicated asymptotic behavior we need to introduce the notion of aj-limit set. 
Let </>' be the flow defined by equation (|3.ip . The w-limit set of x is uj(x) := {y E M. n \3ti, ia--- — * oo s.t. 
4>(tk, x) — > y as k — » oo}. The w-limit sets have the following properties that we will use later. For more 
details see Robinson [9]. 

(i) If cj)(t, x) = y for some t E R, then cu(x) = uj(y). 

(ii) lu(x) is a closed subset and both positively and negatively invariant (contain complete orbits). 

Next we will prove a version of LaSalle's Theorem. For a more general result see the original work 
of LaSalle 0. 

Theorem 3.1. Let xq be an equilibrium point for \3.1}) and U a small compact neighborhood of Xq. 
Suppose there exists L : U — > R a C 1 differ entiable function with L(x) > 0, L(xq) = and L(x) < 0. 
Let E be the set of all points in U where L{x) — 0. Let M be the largest invariant set in E. Then there 
exists a small neighborhood V C U such that u>{x) C M for all x <EV . 

Proof. The conditions of the theorem ensures the stability of Xq. There is a small compact neighborhood 
U of xq and a smaller neighborhood V C U such that <f)(t, x) G U for any x E V and t > 0. As U is 
closed we also have u>(x) C U. 

Let x E V and y £ uj(x) be arbitrarily chosen. Since L(<j)(t,x)) < we have that L(<fi(t,x)) is a 
nonincreasing function of t. Because L is a positive bounded function on U and <p(t,x) remains for all 
time in U we have lim^oo L(<f>(t, x)) — I, where < I < oo. As y E oj{x) and L is continuous we obtain 
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that L(y) = I. The invariance of oj(x) shows that L(<fi(t,y)) = I and consequently L(<j)(t,y)) = for all 
tel. Hence y 6 E and so oj(x) C E. As ui(x) is invariant implies that uj(x) C M. □ 

The following is the main result of this paper. 

Theorem 3.2. Let (R n ,II, H) be a Hamilton-Poisson system and xo G M™ an equilibrium state for the 
dynamic 

x = n(a:) • Vff(a;). (3.2) 

Suppose that there exists a function <p 6 C 00 (IR fe ,K) ; where k equals the number of functionally inde- 
pendent Casimirs for the Poisson structure II, such that 

(i) 5H v {x e ) = 

(ii) 5 2 H ip (x e ) is positive definite, 
where 

H v (x)=H(x)+L P (C 1 (x),...,C k (x)) 

with C\, . . . Ck € C°°(]R n ,IR) a set of functionally independent Casimirs ofH. 

Let G be the matrix defined by (|2.2| then there exist a small closed and bounded neighborhood U of 
the equilibrium state x e of the corresponding metriplectic system 

x = IL(x) ■ VH{x) + G{x) ■ Vtp{C u C k ){x) (3.3) 

and a neighborhood V C U such that every solution of (|3.3|) starting in V approaches U R E as t — > oo, 
where E := {x S W l \VH(x) and V</?(Ci, . . . , Ck){x) are linearly dependent}. 

Proof. First we have to prove that x e is an equilibrium point for the dynamics (13. 3p . This is guarantied 
by Proposition 12.21 

Next we will prove that under the hypothesis of the theorem the function L S C°° (R" , R) given by 

L(x) = Htp(x) - H tp (x e ) 

is a Lyapunov function for the dynamic (|3.3|) . More exactly using the hypothesis and Remark 12. 11 we 
obtain that there exists a closed and bounded neighborhood U of the critical point x e such that 

(i) L{x e ) = 0. 

(ii) L(x) > 0, (V) xeU,x^x e 
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(iii) L(x) < 0, (V) x e U, 

which implies that Xg IS cL stable equilibrium point for (|3.3|) . 

By Remark 1 2. II we have that E equals the set of all points where L(x) — 0. Using Corollary 12.31 we 
have that the largest invariant subset in E for (|3 . 3|) equals E. 

We showed that all the conditions of the Theorem 13.11 are satisfied and so we obtain the desired 
result. 

□ 

Remark 3.1. The above result tells us in fact how to built in an effective way a set controls which locally 
asymptotically stabilize a nonlinear stable equilibrium state of a given Hamilton-Poisson dynamics. 

□ 



4 Examples 



It is well known that Euler's angular momentum equations of the free rigid body can be written on 
in the following form: 

x\ = aix 2 x 3 
x 2 = a 2 xix 3 
x 3 = a- i x 1 x 2 



where 



1 1 



1 1 



1 1 



oi - — - — ; a 2 - — -— ; a 3 - — - — ; 

^3 H 1 l *3 J 2 J-l 



(4.1) 



Iiih, I3 being the components of the inertia tensor and we suppose as usually that 

h > h > h > 0. 

The equations (|4.1[) have the following Hamilton-Poisson realization: 

((*>(3))* «K 3 , {-,.}_, H) 



where {•, •}_ is minus-Lie-Poisson structure on ((so (3))* ~ M. 3 , generated by the matrix: 



n_ = 






- %3 


X 2 


X3 





- X\ 


-x 2 


Xi 
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and the Hamiltonian H is given by: 



. 1 ( x\ x\ Xo 
H(XI,X 2 ,X 3 ) = - + 



2 \h h h 
It is not hard to see that the function C £ C°° (R 3 , R) given by: 

C(xi, x 2 ,x 3 ) = ^(^l + x 2 + £3) 



(4.2) 



(4.3) 



is a Casimir of our configuration ((so (3))* « R 3 , {•, •}_) . 

Let (Mo, 0,0) be an equilibrium point for (|4.1[) . The function H v (x) — H(x) + (p(C(x)), where (p is 
given for instance by 

*(-) = (-i^) a -£, 

satisfies the conditions (i) and («) of Theorem 13.21 In this case the perturbed system (|3 . 3[) is given by 



xi = aix 2 x 3 + X\(x\ + x 2 + X3 — c)( — — x 2 + -7-X3) 

^3 

/ 9 . 9 . 9 \ / 2 ^1 9 \ 

x 2 = a 2 xix 3 + x 2 (xf + x 2 + X3 - c)(— x\ — pXg) 



/1 



/ 3 



(4.4) 



X3 = a 3 xix 2 + x 3 (xf + x 2 + x\ — c)( — ^-x\ + y~xI) 



where 



MS -jr. 



The set of points in R 3 where Vif(x) and V(p(C(x)) are linearly dependent is given by 
E = {(A, 0, 0)|A e R} |J{(0, A, 0)|A e R} (J{(0, 0, A)|A e R}. 



By Theorem l3.2l we obtain that there exists a small closed and bounded neighborhood U C R around 
the equilibrium point (Mo, 0, 0) and V C U such that any solution of (|4.4p starting in V approaches the 
set {(M + A,0,0)|A e [-e,e] C R} as t -> 00. 

A similar result with obvious modifications can be also obtain for the equilibrium state (0,0, Mo), 
M 6 R. 
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